Breaking of factorization of two-particle correlations in hydrodynamics 
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The system formed in ultrarelativistic heavy-ion collisions behaves as a nearly-perfect fluid. This 
collective behavior is probed experimentally by two-particle azimuthal correlations, which are typi- 
cally averaged over the properties of one particle in each pair. In this Letter, we argue that much 
additional information is contained in the detailed structure of the correlation. In particular, the 
correlation matrix exhibits an approximate factorization in transverse momentum, which is taken as 
a strong evidence for the hydrodynamic picture, while deviations from the factorized form are taken 
as a signal of intrinsic, "nonflow" correlations. We show that hydrodynamics in fact predicts factor- 
ization breaking as a natural consequence of initial state fluctuations and averaging over events. We 
derive the general inequality relations that hold if flow dominates, and which are saturated if the 
matrix factorizes. For transverse momenta up to 5 GeV, these inequalities are satisfied in data, but 
not saturated. We find at least as large factorization breaking in event-by-event ideal hydrodynamic 
calculations as in data, and argue that this phenomenon opens a new window on the study of initial 
fluctuations. 



I. INTRODUCTION 

In relativistic heavy-ion collision experiments a large 
second Fourier harmonic is observed in two-particle cor- 
relations as a function of relative azimuthal angle [H-[l| ■ 
This has long been considered a sign of significant col- 
lective behavior [5j , or "elliptic flow" , indicating the ex- 
istence of a strongly- interacting, low-viscosity fluid Q. 
However, only recently has it been realized that all such 
correlations observed between particles separated by a 
large relative pseudorapidity could be explained by this 
collective behavior .7, 8], at least for the bulk of the sys- 
tem. 

One significant piece of evidence for this view was 
the recent observation of the factorization [9j-[l2j of two- 
particle correlations into a product of a function of prop- 
erties of only one of the particles times a function of the 
properties of the second. Specifically, for pairs of parti- 
cles in various bins of transverse momentum Pt-, factor- 
ization of each Fourier harmonic was tested as [9( : 

V nA (Pt,Pt) = {cosn{0 a - 4> b )) I v n {p a T ) x v n (p b T ), (1) 

where the brackets indicate an average over pairs of par- 
ticles (a and b) coming from the same event as well as an 
average over a set of collision events, and 4> a ((j) b ) is the 
azimuthal angle of particle a(b). The left-hand side is a 
(symmetric) function of two variables, p 1 ^ and p\, and in 
general may not factorize into a product of a function v n 
of each variable individually. The fact that this factoriza- 
tion holds at least approximately, then, is a non-trivial 
observation about the structure of the correlation. 

While most known sources of non-flow correlations do 
not factorize at low px [13| , a type of factorization comes 
naturally in a pure hydrodynamic picture where particles 
are emitted independently. They thus have no intrinsic 
correlations with other particles, carrying only informa- 
tion about their orientation with respect to the system 



as a whole. This causes the two-particle probability dis- 
tribution in a single collision event to factorize |14j into 
a product of one-particle distributions, 
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Inspired by this fact, it has often been stated [12|, [l5|, 
[lH that the factorization test in Eq. ((T|) should work 
perfectly in hydrodynamics. The observed approximate 
factorization was hailed as a success for the flow inter- 
pretation of correlations, while small deviations from the 
factorized form was interpreted as a gradual breakdown 
of the hydrodynamic description with increasing trans- 
verse momentum, and of increasing contribution from 
other sources of correlations. 

In this work, we show that factorization as in Eq. ([T]) 
is not necessarily present even in an ideal hydrodynamic 
system governed by Eq. ([2]) . We show that the correlation 
matrix satisfies general inequalities, which are saturated 
by Eq. (QJ. We test these inequalities on ALICE data 
and point out where breaking of factorization occurs. We 
then illustrate with a full event-by-event hydrodynamic 
calculation that the same deviation seen in experiment is 
also present in ideal hydrodynamics. 



II. HYDRODYNAMICS AND TWO-PARTICLE 
CORRELATIONS 



We begin by recalling the discussion originally found 
in Ref. [17| . In a pure hydrodynamic picture, particles 
are emitted independently from the fluid at the end of 
the system evolution according to some underlying one- 
particle probability distribution. One can write any such 
distribution as a Fourier series in the azimuthal angle <p 
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of the particles 
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where = {e m< ^} is the nth complex Fourier flow co- 
efficient, and curly brackets indicate an average over the 
probability density in a single event. Writing V n = 
i> n e m *" , where v n is the (real) anisotropic flow coefficient 
and the corresponding phase, and using V- n = V* 
(where V* is the complex conjugate of V n ), this can be 
rewritten as 

Itt dN \—~ \ 

~N~d6 = 1 + ^ v n (pr,r)) cos n{(f> - ^ n {p T ,V)) ■ (4) 

^ 71=1 

Note that, for this form to describe an arbitrary distri- 
bution, both v n and ^ n may depend on transverse mo- 
mentum px and pseudorapidity T]. 

In this picture, the relation in Eq. ([2]) holds, and a 
complex Fourier harmonic of the two-particle correlation 
factorizes in each event as: 



VZV*=vZv*e«**i-<\ (5) 

The experimental quantity, Eq. ([1]), is then obtained by 
averaging over events: 

V nA (p a T,p b T) = (VnVn*) = (<l£e**<*»-*»>) (6) 

Due to parity symmetry, only the real part remains after 
this average, hence the cosine in Eq. ([TJ. 

From this relation alone, one can make the following 
general statements about the event-averaged correlation 
matrix when flow is dominant: the diagonal elements 
must be positive, and the off-diagonal elements must sat- 
isfy a Cauchy-Schwarz inequality, 

V nA (p a T,P a T) > 0, (7) 

V nA (PT,PT) 2 < V nA (PT,PT)VnA(PT,P b T )- (8) 

Factorization, Eq. (JTJ, implies that the second inequality 
is saturated, i.e., equality is achieved. While flow does 
not necessarily imply factorization, any violation of these 
inequalities is an unambiguous indication of the presence 
of non-flow correlations. 

An inspection of published data from the ALICE Col- 
laboration Q shows that these inequalities are indeed 
violated in certain regimes [181 ] . For n = 3, diagonal ele- 
ments V3a(Pt>Pt) are negative above 5 GeV for 0-10% 
centrality, and above 4 GeV for 40-50% centrality. This 
is a clear indication that there are nonflow correlations 
at high px- For instance, the correlation between back- 
to-back jets typically yields a relative angle A0 ~ n, thus 
producing a negative V3A at high pr- For n = 1, diagonal 
elements are negative not only at high px (with a slightly 
higher threshold than for n = 3), but also for px between 



1 and 1.5 GeV. This is believed to be caused by the corre- 
lation from global momentum conservation (l2l [l9j , but 
it is interesting to note that its effect can be noticed by 
a simple inspection of elements. 

In order to check the validity of the second inequality 
(151), we introduce the ratio 



V nA (p T ,P T ) 



yJv nA (jp T ,p T )V nA {p T ,p b T ) 



(9) 



which is defined when diagonal elements V nA (p T ,p T ) 
and VtcA^T'Pt) are both positive, and lies between — 1 
and +1 if Eq. ([HJ holds. Factorization corresponds to 
the limit r n = ±1. Figure [T] displays T2 and r% as a 
function of p^ and p h T for Pb-Pb collisions at 2.76 TeV, 
0-10% centrality. ALICE results for r-i satisfy the in- 
equalities (jSj) at all pr- When both particles are below 
1.5 GeV, the inequality is saturated, r 2 = 1, within er- 
rors. As soon as one of the particles is above 1.5 GeV, 
however, r 2 is smaller than unity, and the difference with 
unity increases with the difference p^ —p\. Results for r$ 
are qualitatively similar below 5 GeV, with larger error 
bars. However, is closer to 1 than r 2 between 2 and 
3 GeV. The values of r n for mid-central collisions (40- 
50% centrality, not shown) are comparable to the values 
for central collisions, although r 2 is slightly closer to 1. 

The ALICE collaboration concluded from their analy- 
sis that factorization holds approximately for n > 1 and 
Pt below 4 GeV. However, their results actually show 
evidence for a slight breaking of factorization for n = 2, 
as soon as one of the particles has px > 1.5 GeV. Even 
though factorization is broken, the general inequalities 
implied by flow arc satisfied for n = 2 and n = 3 below 
5 GeV for central collisions. It is therefore worth inves- 
tigating in more detail to what extent the breaking of 
factorization which is seen experimentally can be under- 
stood within hydrodynamics. 

First, we recall under which conditions factorization 
holds in hydrodynamics. It implies that the Cauchy- 
Schwarz inequality (jHJ is saturated. By inspection of 
Eq. j6j), this in turn implies that the complex flow vec- 
tors y n a and y„ b are linearly dependent. This is true only 
under the following assumptions: 

1. By parity symmetry, ^ n — ^> b n = in each event. 
I.e., ^ n does not depend on pr, which removes the 
exponential from the right-hand side of Eq. ([5])- 

2. v u {pt) changes from event to event by only a global 
factor, with no pr-dependent fluctuations. v n (pr) 
in the right-hand side of Eq. (QJ then represents the 
rms value over events. 

In general, fluctuations ensure that these conditions 
are not met exactly, and the factorization of Eq. (QJ will 
not be perfect. Within hydrodynamics, the ratio r n in 
Eq. (J9j) has a simple interpretation. Inserting Eq. (|6|) into 
Eq. one obtains 
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FIG. 1. (Color online) Ratio of nondiagonal to diagonal correlations, defined by Eq. versus p?p andpy. Filled stars: ALICE 
data for 0-10% central Pb-Pb collisions at 2.76 TeV [IJ. Open cicles: ideal hydrodynamic calculations for 0-10% central Au-Au 
collisions at 200 GeV. 
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The ratio r n thus represents the linear correlation be- 
tween the complex flow vectors at momenta p?p and p^. 
Since in each event, V£ is a smooth function of pj,, one 
expects that the correlation is stronger when p^ ~ p\, 
and decreases as the difference between p^ and p h T in- 
creases, as a result of the decoherence induced by initial 
fluctuations. ALICE data confirm this qualitative expec- 
tation. 

Note that even in a single hydrodynamic event, factor- 
ization holds in the complex form Eq. (O, but is broken 
if one takes the real part before averaging over particle 
pairs, as in Eq. (JXJ) . The ratio r n in Eq. (fTU|) is then 
cosn(\I/£, — ^n)j which is smaller than unity as soon as 
the flow angle ^ n depends on px- 

The question then becomes: how large arc 
factorization-breaking effects in hydrodynamics, and do 
they have the same properties as seen in data? If purely 
hydrodynamic calculations give the same result as exper- 
iment, then the observed breaking of factorization may 
not indicate the presence of non-flow correlations. 



III. IDEAL HYDRODYNAMIC CALCULATIONS 

To illustrate these concepts we perform calculations 
using the NeXSpheRIO model [20(. This model solves 



the equations of relativistic ideal hydrodynamics with 
fluctuating initial conditions given by the NeXuS event 
generator [2l|. It has proven succesful in reproducing 
RHIC results, in particular the structure of two-particle 
angular correlations in Au-Au collisions at the top RHIC 
energy 22] . It has recently been shown to reproduce the 
whole set of measured anisotropic flow data |23l - [25| . Our 
calculations are therefore performed for Au-Au collisions 
at the top RHIC energy, not for Pb-Pb collisions at LHC 
energy, as would be appropriate for a direct quantitative 
comparison with ALICE data. Our results are merely 
meant as a proof of concept, and as a prediction for mea- 
surements at RHIC. 

We run 3000 hydro events in each 10% centrality 
bin. Anisotropic flow is calculated accurately in every 
event [26| . The ratio r n is displayed in Fig. [1] for n = 2 
and n = 3. Deviations from the factorization limit r = 1 
are already seen at low momentum but become larger 
as the difference between pj, and p\ increases, as ex- 
pected from the general arguments above. Surprisingly, 
the breaking of factorization appears larger in hydrody- 
namics than in experiment. 

The ALICE collaboration has studied factorization 
by performing a global fit of the measured correlation 
Kia(Pt"Pt) k v the right-hand side of Eq. ((T|), where 
v n {pT) is a fit parameter [9]. The ratio of the measured 
correlation to the best fit differs from unity if factoriza- 
tion is broken. We can apply the same procedure to 
our hydrodynamic results. The result is shown in Fig. [3] 
Again, hydrodynamic calculations and experimental data 
show similar trends, with the noticeable difference that 
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FIG. 2. (Color online) Ratio of the left-hand side to the right-hand side of Eq. {T}. Filled stars: ALICE data for 0-10% central 
Pb-Pb collisions at 2.76 TeV [9(. Open cicles: ideal hydrodynamic calculations for 0-10% central Au-Au collisions at 200 GeV. 



the breaking of factorization is significantly stronger in 
ideal hydrodynamics than in data. 

Several effects can explain this discrepancy. First, 
the av erag e px is significantly larger at LHC than at 
RHIC 27], so that it might be more natural to com- 
pare, e.g., 4 GeV at RHIC to 5 GeV at LHC, rather than 
doing the comparison at the same px ■ The second effect 
is viscosity, which is neglected in our calculation. Shear 
viscosity, in particular, tends to damp the effect of ini- 
tial fluctuations |28[. It is therefore natural that it will 



also decrease the breaking of factorization induced by ini- 
tial fluctuations. A similar observation is that the linear 
correlation between the initial eccentricity and the final 
anisotropic flow is stronger in viscous hydrodynamics [29| 
than in ideal hydrodynamics p6| . 



IV. CONCLUSIONS 

We have demonstrated that the detailed structure of 
two-particle angular correlations contains much more in- 
formation than traditional analyses of anisotropic flow, 
where the correlation is averaged over one of the parti- 
cles [3(| ■ Even though such two-dimensional analyses are 
much more demanding in terms of statistics than tradi- 
tional analyses, they bring new, independent insight into 
the underlying physics. 

In particular, we have shown that quantum fluctua- 
tions in the wavefunction of incoming nuclei result in 
a decoherence in the angular correlations produced by 
collective flow, which becomes increasingly important as 



the difference between particle momenta increases. Due 
to this effect, factorization of angular correlations is bro- 
ken even if collective flow is the only source of correla- 
tions. Our numerical calculations show that factoriza- 
tion breaking is at least as strong in ideal hydrodynam- 
ics as in experimental data, thereby suggesting that all 
correlations below pt ~ 5 GeV (for central Pb-Pb colli- 
sions at LHC near midrapidity) are actually dominated 
by flow. The sensitivity of this decoherence phenomenon 
to viscosity has not yet been investigated, but we antic- 
ipate that factorization should be restored as viscosity 
increases, thus potentially offering a new means of con- 
straining the viscosity from data. 

Decoherence also provides a natural explanation for 
the important observation that event-by-event fluctua- 
tions reduce elliptic flow at high px [31], thus improv- 
ing agreement between hydrodynamics and experimen- 
tal data. Indeed, v-i at high pt is inferred from az- 
imuthal correlations between a high px particle and all 
other particles — mostly low px particles, and these az- 
imuthal correlations are reduced due to the decoherence 
phenomenon. Note that the other main explanation for 
the reduction of vi at high px, viscosity, typically relies 
on the assumption of a quadratic momentum dependence 
of the viscous correction to the distribution function at 
freeze-out Sf, which may not be correct [32| . 

In this paper, we have focused on the transverse mo- 
mentum dependence of the correlations. The rapidity 
dependence of the correlation is also worth investigating. 
In particular, it was recently observed that azimuthal 
correlations decrease as a function of the relative pseu- 
dorapidity [33] , at variance with common lore that corre- 
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lations due to flow are essentially independent of rapidity. 
While standard models of initial conditions do predict a 
mild rapidity dependence of azimuthal correlations [13], 
longitudinal fluctuations [35| could also produce a deco- 
herence effect similar to the one studied here. The de- 
tailed structure of two-particle correlations as a function 
of both particle momenta thus opens a new window on 
the study of flow fluctuations. 
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